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We present observations of electromagnetic precursor 
waves, identified as whistler mode waves, at supercritical 
interplanetary shocks using the Wind search coil magne- 
tometer. The precursors propagate obliquely with respect 
to the local magnetic field, shock normal vector, solar wind 
velocity, and they are not phase standing structures. All 
are right-hand polarized with respect to the magnetic field 
(spacecraft frame), and all but one are right-hand polar- 
ized with respect to the shock normal vector in the normal 
incidence frame. They have rest frame frequencies f ci < f 
<C fee and wave numbers 0.02 < k p ce < 5.0. Particle dis- 
tributions show signatures of specularly reflected gyrating 
ions, which may be a source of free energy for the observed 
modes. In one event, we simultaneously observe perpendic- 
ular ion heating and parallel electron acceleration, consis- 
tent with wave heating/acceleration due to these waves. Al- 
though the precursors can have <SB/B 0 as large as 2, fluxgate 
magnetometer measurements show relatively laminar shock 
transitions in three of the four events. [Date: 03-27-2012] 


1. Introduction 

Collisionless shock waves, predicted over 50 years ago 
[ Petschek , 1958], are a ubiquitous nonlinear structure ob- 
served in nearly every space plasma environment. Shock 
waves require energy dissipation - transfer of energy from 
one form to another through an irreversible process - where 
the type of energy dissipation in collisionless shock waves 
depends upon the Mach number, My, and has implica- 
tions for the global structure of the shock wave [Coroniti, 
1970]. Theory and observations suggest that low Mach num- 
ber shocks rely upon wave dispersion [Kennel et al., 1985] 
and/or anomalous resistivity due to wave-particle interac- 
tions for energy dissipation [ Sagdeev , 1966; Wilson III et al, 
2007]. Above the first critical Mach number, M cr , the shock 
requires additional energy dissipation in the form of particle 
reflection to limit wave steepening [Edmiston and Kennel, 
1984], 

Observations of low Mach number quasi-perpendicular 
(shock normal angle, d Bn , > 45°) low-/! shocks show a rel- 
atively laminar transition, with a compressive electromag- 
netic (EM) precursor wave supplying the dispersive dissi- 
pation [Kennel et ai, 1985]. The precursor wave has been 
identified as a right-hand polarized wave with f ci < f -C ice, 
where f cs is the cyclotron frequency of species s and i(e) 
represent ions(electrons), consistent with a whistler mode 
wave. These ’whistlers’ have been observed in nearly ev- 
ery space plasma environment ([ e.g . Wilson III et al., 2009, 
and references therein]), can couple to multiple wave modes 
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[Marsch and Chang, 1983; Wu et al, 1983; Matsukiyo and 
Scholer, 2006], and can interact with both ions and electrons 
[Wu et al., 1983]. Recent PIC simulations have found that 
whistler precursors may play an important role in particle 
acceleration at higher Mach number shocks [e.g. Riquelme 
and Spitkovsky, 2011]. Thus, whistlers are thought to play 
an important role in shock wave dynamics through both dis- 
persion and anomalous resistivity. 

In this study, we present the first search coif observa- 
tions of whistler precursors in and around the ramp region 
of four supercritical interplanetary (IP) shocks observed by 
the Wind spacecraft. We also present evidence for wave 
heating and particle acceleration in one event. This work 
presents the first 3-D magnetic field waveform observations 
of waves with f sc < i lh and i lh < f sc < f ce at supercritical IP 
shocks, where f sc is the spacecraft frame frequency and Ilh 
(= VYeeTci) is the lower hybrid resonance frequency. 

2. Observations 

Waveform captures were obtained from the Wind /WAVES 
instrument [Bougeret et al, 1995], using the time domain 
sampler slow (TDSS) receiver, providing three magnetic 
field components and one electric field component, with 2048 
data points sampled at 1875 samples/s. The ambient mag- 
netic field was obtained from the dual, triaxial fluxgate mag- 
netometers [Lepping et al., 1995] sampled at ~11 samples/s. 
Ion and electron moments obtained from the Wind/3DP 
EESA and PESA particle detectors [Lin et al, 1995] were 
used in conjunction with shock parameters found at [Kasper, 
2007] to calculate M cr , assuming a polytrope index, 7 = 5/3 
[Edmiston and Kennel, 1984]. Supplemental ion distribu- 
tions were calculated from the two Faraday Cup (FC) ion 
instruments from the Wind/S WE experiment [Ogilvie et al, 
1995]. For more details about the WAVES and 3DP instru- 
ments and analysis, see Wilson III et al. [2010]. 

The top half of Table 1 lists the relevant shock param- 
eters, provided by [Kasper, 2007], for the four events pre- 
sented herein, which include: the shock normal speed in the 
spacecraft frame, V s(ln ; the angle between the shock normal 
vector and upstream magnetic field, 9 Bn \ fast mode Mach 
number, My,; and the shock compression ratio N i2 /N tl . The 
bottom half shows the ratios of the My to M cr as well as 
to the three whistler critical Mach numbers [Krasnoselskikh 
et al, 2002]: M„ corresponds to the maximum Mach num- 
ber at which a linear whistler can phase stand with respect 
to the shock; M gr is the maximum Mach number at which 
the wave can carry energy into the upstream; and Wl nw is 
the maximum Mach number for which a stationary solution 
can be found above which the wave breaks. All four events 
have My/M cr > 1 (i.e., are supercritical) and two of the 
events have My/M nra > 1. 

3. Analysis 

Figure 1 plots the Wind/MFI data at ~11 samples/s for 
the four supercritical IP shocks examined. The only event 
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with resolvable magnetic fluctuations is shown in Figure 1A. 
The other three shock profiles show relatively laminar tran- 
sitions, and magnetic foots are noticeable in Figure IB and 
ID. All four IP shocks exhibit magnetic overshoots, consis- 
tent with previous observations of supercritical shocks [ Far- 
ris et al. , 1993]. Note that Figure 1A was examined by 
Wilson III et al. [2009], but only with fluxgate magnetome- 
ters. The important observation here is that the magnetic 
fields are under-sampled and that higher frequency fluctu- 
ations (discussed and shown below) with amplitudes com- 
parable to the shock ramp cannot be observed in fluxgate 
magnetometer data at this sample rate. High time resolu- 
tion measurements (>40 samples/s) are needed to resolve 
the ramp as one can see the spiky nature of these plots is 
indicative of undersampling. 

Figure 2 shows the four TDSS samples that span the 
ramps of the four supercritical IP shocks. Each panel shows 
the three normal incidence frame (NIF) components [defined 
by Sundkvist et al., 2012] of the search coil magnetic field. 
The ratio of the maximum peak-to-peak wave amplitude to 
the ambient magnetic field at the start of each event, <5B/B o; 
is: (A) ~0.9; (B) ~2.1; (C) ~2.0; and (D) ~0.5. Thus, at 
least for B and C, the wave amplitudes significantly modify 
the ramp structure shown in Figure 1. 

The waveforms were analyzed using Minimum Variance 
Analysis (MVA) [ Khrabrov and Sonnerup, 1998] to deter- 
mine the wave vector, k; polarization with respect to the 
ambient magnetic field; and wave normal angles with re- 
spect to the local magnetic field (0 kB )\ shock normal vector 
(#*,„); and local solar wind velocity (9 kv ). Wavelet analysis 
revealed that these waves exhibit fluctuations at both high 
(fih < < fee) and low (f sc < f lh ) frequencies. Therefore, 

we applied multiple standard Fourier bandpass filters to each 
waveform prior to using MVA on specific subintervals (see 
Wilson III et al. [2009] for more details), requiring that the 
ratio of the intermediate-to-minimum eigenvalues be >10 to 
use a result. This analysis resulted in 84 unique k-values. 

All four precursors are right-hand polarized with respect 
to the ambient magnetic field, B 0 , in the spacecraft frame of 
reference. Only the precursor shown in Figure 2B shows a 
left-hand polarization in the NIF with respect to the shock 
normal vector, n, the rest are right-handed. The precursors 
exhibit broad range of values for 6 kB , 6 kn , and 9 kv , con- 
sistent with theory [Wu et al., 1983] and prior observations 
[Wilson III et al., 2009; Sundkvist et al., 2012]. For the low 
frequency components, the majority had 30° < 9 kB < 75°, 
15° < 9 kv < 45°, and 15° < 9 kn < 60°, consistent with pre- 
vious observations of non-phase standing precursors at IP 
shocks [Wilson III et al., 2009]. For the high frequency com- 
ponents, the majority had 9 kv and 9 kn > 45°, while 9 kB was 
broad from ~10° to 90° with most < 60°. Their right-hand 
polarization with respect to B„, propagation directions, and 
frequencies are all consistent with oblique whistler modes in 
the solar wind. 

Since we know k and the solar wind velocity vector is well 
defined at the start of each TDSS sample, we can use the 
Doppler-shifted cold plasma dispersion relation for whistlers 
to estimate the magnitude of the wave vector [see Coroniti 
et al., 1982]. We found the waves to have 0.03 < kc /oo pe < 
4.0 or 0.02 < k p ce < 5.0, where c is the speed of light, u> pe is 
the electron plasma frequency, and p ce = V Te /w pe is the elec- 
tron thermal gyroradius (where V Te = \Z2k B T e /m e is the 
electron thermal speed, T e the electron temperature, and k B 
is Boltzmann’s constant). Using these results, we estimated 
that the precursors had rest frame frequencies 0.01 < f/f lh 
<10. These results are consistent with previous solar wind 
observations [Coroniti et al, 1982] and theory /simulations 
[Marsch and Chang, 1983; Wu et al., 1983; Lembege et al., 


2009] of whistlers. At these frequencies, whistlers lie on the 
same branch of the dispersion relation as both EM lower 
hybrid [Marsch and Chang, 1983] and magnetosonic wave 
modes [Wu et al., 1983]. 

We examined the particle distribution functions (Figures 
3i and 3e) from the Wind/3DP instrument for evidence of 
free energy or signatures of ion reflection. Only three of the 
events (see Figures 1A, IB, and ID) had the burst parti- 
cle data (full distribution every 3 seconds) needed to exam- 
ine the distributions in and around the shock ramps. We 
were able to produce the electron distributions using ~11 
sample/s MFI data, reducing the aliasing across the shock, 
similar to method introduced by Schwartz et al. [2011]. This 
was not possible for the PESA distributions. However, when 
we compared the distributions produced using multiple in- 
stantaneous magnetic field vectors to those produced using 
only one (typical method), the results were consistent with 
each other for the Figure IB event. This is primarily due 
to the field direction being predominantly in the same di- 
rection throughout the duration of the distribution. Thus, 
we believe that the PESA High distributions produced using 
only one magnetic field vector (Figure 3i) accurately repre- 
sent the ion distributions across the shock in Figure IB. 

All three events with burst particle data showed evidence 
for gyrating ions (example in Figure 3i) at/near the pre- 
dicted gyrospeed for specular reflection [Thomsen et al., 
1983a]. The particle distributions in Figure 3 have color- 
coded outlines that correspond to the color-coded regions 
in Figure IB. The IP shock associated with Figure IB also 
showed evidence for strong perpendicular (with respect to 
the magnetic field) ion heating (second panel of Figure 3i), 
which was supported by PESA Low and SWE observations. 
The ion heating is apparent from ~170-700 km/s (or ~150- 
3000 eV proton) . This anisotropic heating was not observed 
in the two other events with burst particle data, yet the 
magnetic compression ratio is similar for each event (B 2 /Bi 
« 2. 1-2.7, see Figure 1). Therefore, we argue that the ob- 
served differences in ion heating are due to the larger waves 
observed in Figure 2B. The observed gyrating ions may pro- 
vide the free energy for the modes observed in Figure 2, con- 
sistent with theory [Wu et al., 1983]. 

Figure 3e shows the three electron distributions, pro- 
duced using ~11 sample/s MFI data, at the same times 
as the ion distributions in Figure 3i. The electron distri- 
butions show features consistent with previous observations 
near collisionless shocks [e.g. Thomsen et al., 1983b]. The 
middle panel shows two interesting features as evidence for 
parallel electron acceleration. The first, at low energy, is 
the “bump-on-flattop” feature that is likely due to accel- 
eration of the thermal core by the quasi-static cross shock 
potential [Thomsen et al., 1983b]. The second, at higher 
energy, is a shoulder-like feature (shaded magenta region) 
that has been predicted in multiple simulation studies as 
the signature of highly oblique waves stochastically accel- 
erating electrons [e.g. Wu et al, 1983; Cairns and McMil- 
lan, 2005]. The same studies show that these waves also 
cause perpendicular ion heating of the lower energy ions, 
consistent with our observations. We also observe higher 
frequency electrostatic waves near these events, which may 
explain the flattop distribution in the third panel of Fig- 
ure 3e. However, these waves, which have been examined 
previously [Wilson III et al., 2007], cannot explain both the 
perpendicular ion heating and parallel electron acceleration. 
Therefore, we argue that the wave observed in Figure 2B is 
responsible for the strong anisotropic ion heating in Figure 
3i and parallel electron acceleration in Figure 3e, consistent 
with the predicted interaction between low frequency (i.e. 
<f I(l ) very oblique (i.e. 9 kB > 75°) whistler waves and par- 
ticles [Wu et al., 1983]. 
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Note that the shoulder-like feature in Figure 3e is di- 
rected toward the downstream. This coupled with the wave 
in Figure 2B having a left-hand polarization in the NIF 
with respect to n are consistent with a downstream directed 
phase velocity, consistent with recent instability driven sim- 
ulation results [Comi§el et al, 2011]. The rest of the waves 
in Figure 2 were consistent with recent observations [Sund- 
kvist et al, 2012] that suggest the precursors are dispersively 
driven. Though the wave in Figure 2B appears to have a 
downstream directed phase velocity, it may still have an up- 
stream directed group velocity. For instance, the dispersion 
relation for oblique magnetosonic waves [e.g. Coroniti , 1970] 
gives a range of wave numbers satisfying 0.5 < (kc)/cu pe < 
1.0 where the wave in Figure 2B could have upstream di- 
rected group velocities but downstream directed phase ve- 
locities, which would be consistent with the results of Sund- 
kvist et al. [2012]. The downstream directed phase velocity 
suggests that the wave in Figure 2B is receiving some free 
energy from reflected gyrating ions that are seen in Figure 
3i. 

4. Conclusions 

We present the first observations of whistler precur- 
sor waves at IP shocks above both the first and nonlin- 
ear whistler critical Mach numbers. The precursors have 
large amplitudes (<5B/B„ > 2), have right-hand polariza- 
tions (spacecraft frame) with respect to B„, and propagate 
obliquely with respect to B„ and n; thus they are not phase 
standing. Only one event showed a left-hand polarization 
with respect to the n in the NIF. The waves have rest frame 
frequencies f ci < f -C f ce and wave numbers 0.02 < k p ce < 
5.0. Their characteristics are consistent with EM whistlers 
and EM lower hybrid waves, both modes known to interact 
with both ions and electrons [Wu et al., 1983]. 

In one event we simultaneously observe with the wave 
specularly reflected gyrating ions and electrons showing sig- 
natures of parallel acceleration directed downstream. This 
event was the only precursor that had a left-hand polariza- 
tion (in NIF) with respect to n. These two features are 
consistent with downstream directed phase velocity, consis- 
tent with recent simulation results [Comigel et al, 2011]. 
Therefore, we argue that the strong perpendicular ion heat- 
ing and parallel electron acceleration observed was due to 
interactions of the particles with this mode. 

Two of the IP shocks in our study, the 1997-12-10 (A) and 
1998-09-24 (C) events, satisfy M ; /M„ < 1 and M. f /M. nw < 
1. Thus, observation of whistler precursors at these two 
IP shocks is consistent with the theory that phase standing 
waves limit the dispersion of the shock [Krasnoselskikh et al, 
2002] and previous observations [Farris et al., 1993]. How- 
ever, the 1998-08-26 (B) and 2000-02-11 (D) events have 
My/Mnu, > 1, which means these two IP shocks should not 
exhibit whistler precursors and should undergo reformation. 
Non-stationary reformation has been observed at the ter- 
restrial bow shock [Lobzin et al., 2007] and in simulation 
[e.g. Comigel et al., 2011]. Since the waves observed in our 
study are not phase standing, it may account for the obser- 
vations when M^/Mnu, > 1. However, unlike the dispersion 
driven reformation discussed by Krasnoselskikh et al. [2002], 
the downstream directed phase velocity observed herein sug- 
gests mechanisms consistent with instability driven reforma- 
tion observed by Comigel et al. [2011]. We note that the ap- 
parent laminar shock transitions observed in Figures 1B-D 
may only be a consequence of undersampling, not a physical 
transition to a supercritical structure as previous observa- 
tions suggested [Farris et al., 1993]. We hypothesize that 
given the right instrumentation, a collisionless shock ramp 
would be observed to be buried in large amplitude electro- 
magnetic turbulence extending from nearly DC frequencies 


up through f ce . 

We have presented the first waveform observations of 
waves near and above f lh simultaneously observed with par- 
ticle distributions showing the expected signatures of inter- 
actions with these wave modes. This study provides further 
evidence that microinstabilities and turbulence can play a 
significant role in particle dynamics and energy dissipation 
at collisionless shocks. 
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Table 1. Shock Parameters and Mach Number Ratios 


Date 

V sh n (km/ s) 

e Bn 

M/ 

N„/N tl 

1997-12-10 

391 ± 12 

71° ± 2° 

2.3 ± 0.1 

2.5 ± 0.4 

1998-08-26 

687 ± 27 

82° ± 3° 

4.7 ± 0.2 

2.9 ± 0.3 

1998-09-24 

772 ± 96 

82° ± 2° 

2.9 ± 0.1 

2.2 ± 0.4 
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Figure 1. Plot of 12-second windows of the Wind flux- 
gate magnetometer data for the four IP shocks. The top 
panel for each event plots the magnitude of the magnetic 
field and the bottom panel plots the three NIF compo- 
nents of the magnetic field defined by Sundkvist et al. 
[2012]. Each event has three vertical lines: the center 
of the shock ramp (green); and the start (red) and end 
(blue) times of the associated TDSS sample in Figure 2. 
The three color-coded shaded regions in panel B show the 
time ranges of the particle distributions shown in Figures 
3i and 3e. 
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Figure 2. Plot of the four unfiltered TDSS sam- 
ples, identified as magnetosonic-whistler precursor waves, 
overlapping with the four IP shock ramps. Each panel 
plots the three NIF components of the magnetic field. 
The vertical green line in each panel marks the center of 
the corresponding IP shock ramp in Figure 1. The rel- 
ative amplitudes of the waveforms are illustrated by the 
vertical black arrows. 
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Figure 3. The evolution of the particle distributions, 
ions (i) on left and electrons (e) on right, across the shock 
ramp in Figure IB. The contour plots show contours of 
constant phase space density in the plane containing the 
ambient magnetic field (horizontal axis of contours) and 
solar wind velocity. Projected onto this plane are the 
following: shock normal direction (dashed red line, both 
i and e); intersection of the shock plane (solid magenta 
line, only e); solar wind velocity (solid black line, both 
i and e); sun direction (solid blue line, only i); heat 
flux vector (solid blue line, only e); and a circle (dashed 
black, only i) at ~700 km/s corresponding to the pre- 
dicted gyrospeed of specularly reflected ions [ Thomsen 
et ai, 1983a]. The right column for each species plots 
the parallel (solid red line), perpendicular (dashed blue 
line), and parallel one-count level (solid green line) cuts 
of the distributions. In the cuts for i, horizontal black 
dotted lines correspond to the circle in the contour plots. 
In the cuts for e, labels defining the respective regions 
and what we argue to be evidence of parallel electron 
acceleration are shown. 


